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We address the tunneling conductance and spin inelastic tunneling spectroscopy of localized para¬ 
magnetic moments in a superconducting environment, pertaining to recent measurements on Fe- 
octaethylporphyrin-chloride using superconducting scanning tunneling microscopy. We demonstrate 
that the Cooper pair correlations in the tip and substrate generate a finite uniaxial anisotropy field 
acting on the local spin moment, and we argue that this field may be a source for the observed changes 
in the conductance spectrum for decreasing distance between the scanning tunneling tip and the local 
magnetic moment. We make a side-by-side comparison between the superconductor-superconductor 
junction and normal-metal-superconductor junction, and find qualitative agreement between the two 
setups while quantitative differences become explicit. When simulating the effects of electron pump¬ 
ing, we obtain additional peaks in the conductance spectrum that can be attributed to excitations 
between higher-energy spin states. The transverse anisotropy field couples basis states of the local 
spin which opens for transitions between spin states that are otherwise forbidden by conservation of 
angular momentum. Finally, we explore the influences of temperature, which tend to enable in-gap 
transitions, and an external magnetic field, which enables deeper studies of the spin excitation spec¬ 
trum. We especially notice the appearance of a low and high excitation peak on each side of the main 
coherence peak as an imprint of transitions between the Zeeman split ground states. 

PACS numbers: 74.55.+v,73.20.Hb,71.70.Gm,75.10.Dg 


I. INTRODUCTION 

Research into single-spin manipulation remains one 
of the most active areas in materials science. This is 
justifiable, as control of single spins would enable infor¬ 
mation storage with an order-of-magnitude increased 
density as well as the possible realization of practical 
quantum computers. 1/2 Writing and reading information 
from a localized atomic or molecular spin necessarily in¬ 
volves controlled transitions between different energy 
states. In most experimental cases under ambient condi¬ 
tions, however, spontaneous interaction with surround¬ 
ing spin carriers severely limits the mean free lifetime of 
the local spin excitations below a realistic clock cycle. 

Magnetic atoms or molecules resting on a metal sur¬ 
face, for example, typically deexcite within picoseconds 
as energy and angular momentum are transferred to itin¬ 
erant electrons of the substrate. 3-5 As a measure to in¬ 
crease such lifetimes, by limiting the number of ways the 
local spin can give away energy, an insulating layer can 
be applied between the metal and the local spin. CuO, 
BN, and CU 2 N have all been used in this manner, to ef¬ 
fectively create a gapped substrate, which increases the 
mean lifetime to hundreds of picoseconds. 3 * 1 s 

While several novel ways to increase spin excitation 
lifetimes have been suggested and proven successful, a 
natural progression from a separating insulating layer is 
to use a superconducting substrate that exhibits a per¬ 
fect band gap yet still conducts charge. 9,10 At low tem¬ 
peratures, a spontaneous deexcitation of the local spin 
state must then provide enough energy to break up a 
Cooper pair in order for the main deexcitation mech¬ 


anism to occur, i.e., quasiparticle-hole pair creation. 6 
A drawback of a superconducting substrate is the ap¬ 
pearance of unwanted Shiba states, within the super¬ 
conducting gap, generated by exchange interaction be¬ 
tween the localized spin moment and the electrons in 
the superconductor. 11-111 To minimize the effect of these 
states Heinrich et al. successfully utilized a paramagnetic 
organic molecule, e.g., M-octaethylporphyrin-chloride 
(M-OEP-C1) where M denotes a transition-metal ele¬ 
ment (Mn, Fe, Co, Ni, Cu), to encage the local mag¬ 
netic moment such that the direct interaction is kept to a 
minimum. 16,1 ' The Shiba states then migrate close to the 
main coherence peaks and are indiscernible unless the 
temperature is very low. In addition to providing sep¬ 
aration, the ligand cage of the paramagnetic molecule 
also generates an environment of magnetocrystalline 
anisotropy for the central magnetic moment that splits 
up the otherwise degenerate spin states into different 
energy levels.' 13 This method prolonged the mean life¬ 
time of the first excitation to t ~ 10ns, which is enough to 
clearly observe pumping into higher spin states. 16 The 
experiments, conducted on a Pb substrate using a Pb 
covered tip at 1.2 K, shows, in addition, that inelastic 
scattering between the tunneling electrons and the local 
spin moment only give signatures in the dl/dV spectra 
for bias potentials |eV| = A su b + A t ; p + A mn , where A mn is 
the spin state excitation energy. 

The theoretical model derived in this paper emulates 
single-electron tunneling in a scanning tunneling micro¬ 
scope (STM) setup where the tip is made up of a normal 
metal (NM) or a superconductor (SC). As a substrate, 
on which a paramagnetic organic molecule lies, only 
a superconductor is considered. The magnetic center 
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FIG. 1: (Color online)Schematic illustration of the STM setup. 
The tip is modeled both as a normal metal and a superconduc¬ 
tor, while only a superconducting substrate is considered. As 
indicated by the arrows, electrons may tunnel directly between 
tip and substrate or through intermediate interaction with the 
local spin S, of the paramagnetic molecule, under bias voltage 
V. 


of the molecule provides a local spin moment, within 
an anisotropic environment, elevated enough to prevent 
significant direct magnetic interaction with the close-by 
superconductors. An applied bias voltage, which con¬ 
trols the relative Fermi levels of the tip and substrate, 
will induce a tunneling current of electrons that either 
pass the local spin moment unnoticed or interact with 
exchange of energy and angular momentum. See Fig. 1 
for a sketched illustration of the setup. 

In excellent agreement with experiment, our transpar¬ 
ent (differential) conductance expression yields signa¬ 
tures of inelastic spin transitions only outside of the tip 
and substrate superconducting gap at low temperatures. 
We also reproduce the observed effects of pumping to 
reveal interactions with higher spin states. Beyond the 
reproduction of experimental results, the conductance 
spectra are thoroughly investigated with respect to vary¬ 
ing anisotropies and external magnetic fields. 

In particular, we investigate the effect of the Cooper 
pair correlations in the tip and substrate on the spectrum 
of the local magnetic moments. We show that these 
give rise to an additional contribution to the uniaxial 
anisotropy, something which can partially explain the 
observed spectral changes upon bringing the STM tip 
closer to the paramagnetic sample. 16,19 

While some of the results regarding the conductance 
spectrum for the SC-SC junction have been published 
elsewhere, 20 the present paper also includes the above- 
mentioned investigation of the Cooper-pair-induced 
uniaxial anisotropy field, the case of a localized magnetic 
moment embedded in a normal-metal-superconductor 
junction, and a systematic study of the influence of tem¬ 
perature and external magnetic fields. For completeness 
and to enable a comparison between the different sce¬ 
narios, we also include details of the SC-SC junction. 

The paper is organized as follows. In Sec. II we define 


the microscopical model for our setup, in Sec. Ill we elu¬ 
cidate the impact of the tip and substrate electrodes and 
the tunneling current on the spin excitation spectrum, in 
Sec. IV we derive the expressions for the tunneling cur¬ 
rent including the expression pertaining to the inelastic 
electron tunneling spectroscopy (IETS) measurements, 
in Sec. V we present and analyze the main results for 
spin S = 1 and S = 5/2 systems, and we conclude the 
paper in Sec. VI. 

II. THEORETICAL DESCRIPTION OF THE MODEL 

The electronic composition and interplay within the 
STM device is governed by the total Flamiltonian 

=< 7"ftip + 'td su b + ( H\ + 'Hs + <H K[ + 'Hks, (1) 

where ; p and 'H su b give the electronic structure of the 
tip and the substrate, respectively. In this study, we 
consider two different types of tip states: (i) normal metal 
and (ii) superconducting. These are modeled using 

, h (NM ] = ^ £p<j£ d r c pCT , (2a) 

PIT 

^tip J = ^ £ v° c po c po + ^ / ^ ti P c pT c -p,i + H' c ' (2b) 

pa pa 

respectively. Here, c pa ( c pa ) creates (destroys) an elec¬ 
tron/quasiparticle in the tip with momentum p and spin 
a =T,|. The substrate quasiparticles are modeled by 

^sub = ^ ^ ^ ^ \ hsubWb 6 _k,[ bl.C., (3) 

ku ka 

where c 1 ^ (cy,,) and k denote electron operators and mo¬ 
mentum, respectively. The three parts, tip, substrate, 
and sample, are connected via tunneling through 

‘Ht = £ [T 0 <W + Tkw ■ S] c kaV + H.c.. (4) 

pk 

ao' 

One direct tunneling path, with a rate Tq (6 aa > is the 
Kronecker delta), is spin preserving, whereas a second 
tunneling path, with rate T\, accounts for the interaction 
between the electron spin and localized spin. Here, <j a o' 
is the Pauli-matrix vector. The ratio T| /Tq may be of 
the order of unity since this rate is determined both by 
the tunneling overlap as well as the Coulomb assisted 
tunneling rate; see, e.g.. Ref. 21. 

In addition to the tunneling contribution “Hj, we con¬ 
sider the effects of Kondo like coupling between the local 
spin moment and the electrons in the tip and substrate, 
respectively. This is introduced through the contribu¬ 
tions 

Kt "t ^Ks ^ I 1 ^ CpO^OO'Cp'ty'a ' Is ^ ^ A CJ',:, y' ; j ‘ S, 

aa’ pp' kk' 

( 5 ) 
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where T t i s \ is the energy for the coupling between the 
local spin moment and the electrons in the tip (substrate). 

The local magnetic moment, embedded in the 
anisotropic environment of the organic molecule, de¬ 
rives its 2S + 1-fold spectrum of spin eigenenergies and 
states {E a ,\a)} from the Hamiltonian 

<H s = -gLi B B-S + DS 2 z + ^(S 2 + +S z _). (6) 

Here, g is the gyromagnetic ratio, lib is the Bohr mag¬ 
neton, and B is an external magnetic field. - For ar¬ 
bitrary integer total spin moment S and a finite uni¬ 
axial anisotropy field D, the basis states \S z ,m z ), rn 7 = 
—S z ,—S z + 1 ,...,S 2 , remain eigenstates with twofold de¬ 
generate excitations. The transverse anisotropy field E 
will split up these excitations as well as cause the eigen¬ 
states to form linear combinations of the basis states. 
A half-integer spin moment behaves much in the same 
way under finite anisotropies, with the exception that £ 
no longer splits up the two fold degenerate excitations 
but rather shifts the energy levels somewhat. 

Direct local interactions with the substrate, respon¬ 
sible for Shiba states, are described by < H^ S . For a 
large local spin moment (S » 1) that couples weakly 
(T s —> 0) to the surface electrons (ensured by the separat¬ 
ing ligand cage), in-gap resonances appear at energies 
a> 0 = ±A sub [l - (nNT s S/2) 2 ]/[l + (nNT s S/2) 2 ], 15 where N 
is the substrate density of states (DOS). The energies 
of the Shiba states, hence, approach the edges, or co¬ 
herence peaks, of the superconducting gap whenever 
(■ nNJS/2) 2 « 1. 

In a quantum mechanical treatment, the energies of 
the Shiba states have to be considered using perturba¬ 
tion theory. To the lowest-order approximation, the en¬ 
ergies of the Shiba states are given by A sub [l - a(N /S/2) 2 ], 
where a = 9 (a = 1) for antiferromagnetic (ferromagnetic) 
coupling /. Despite the different form of the excitation 
energies obtained in this perturbational approach, the 
general statement that the energies approach the edges 
of the superconducting gap in the limit (tcNJS/2) 2 <sc 1 
remains valid. 

We hence conclude that for NJS < 2k bT/ 3n, the Shiba 
are hidden within the thermally broadened coherence 
peaks. Although these states may influence the lifetime 
of the spin excitations, we omit this contribution in the 
following discussion in order to focus on the signatures 
observed in experiment. 

III. MODIFICATIONS TO THE SPIN SPECTRUM 
INDUCED BY TIP, SUBSTRATE, AND CURRENT 

Here we discuss a possible reason and explanation for 
the strong changes in the magnetic anisotropy, acting on 
the local magnetic moment as a function of the tunneling 
current, which was observed in Refs. 16,19. 

From the microscopical model introduced in Eq. (1), 
we can derive an effective spin Hamiltonian which is 


suitable for studies of possible sources of broadening and 
level shift in the spin excitation spectrum. Employing the 
methods used in Refs. 23-26, we arrive at the effective 
model for the local spin, 

( H e s ff =<H S + DS-S + 2 TSl, (7) 

where both fields D and F comprise three contributions 
arising from the local spin interactions with electrons in 
the tip (’Dt/T't), substrate (D s /F s ), and tunneling current 
CDd^Fc)- We notice that the second contribution to the 
spin spectrum, !DS ■ S = DS 2 , does not provide any essen¬ 
tial changes to the relative level spacing in the spectrum, 
except for an overall shift of the energy. As we shall 
see, however, this field contributes a source of dissipa¬ 
tion whenever the voltage bias is greater than 2| A|, which 
generates a current-induced broadening of the spin exci¬ 
tations. Second, we notice that the presence of the super¬ 
conducting electrodes and tunneling current generates 
a contribution to the axial anisotropy in the last term of 
Eq. (7), which was also observed in experiments. 16 ' 19 

The contribution from the current mediated interac¬ 
tions is given by 

£> C =-Tc 

2T 2 y r Gp(cc)G>(a/) - G>(<u)G<(a/) dev dev' 

1 J ev-ev' +eV -ib 2n 2n ' 

pk 

(8a) 

T ,_ 1T iy r n'wnw)-rywtw ^ 

1 ^r J J co-co'-eV-ib 2n 2n 6 

P k 

(8b) 

The tip (substrate) contribution is obtained analogously 
in the limit V —> 0 by replacing k (p) with p' (k'). 
Here, we have introduced the lesser anomalous sur¬ 
face Green functions GFs F<(a>) = i( c kfcr c -kia)( M ) an d 
F^ < (cv) = )(&;), the corresponding greater GFs 

FkM = H)<c—kjCkTaX^) and F^Jcv) = HKc^c^Xcu), 
as well as the corresponding lesser and greater anoma¬ 
lous GFs for the tip in the case of a superconducting 
tip. Otherwise, these propagators are absent for the tip, 
which implies that Ft/ C — 0. 

Before moving ahead, we notice, however, that the 
field F c , which is generated by the Cooper pair corre¬ 
lations in the tunneling current, acquires a harmonic 
temporal variation with the ac Josephson current and, 
hence, vanishes on average. The spin excitation spec¬ 
trum is, therefore, expected to fluctuate around its static 
spectrum with a frequency given by the Josephson cur¬ 
rent. For the IETS measurements in this paper, we will 
therefore neglect the influence of this anisotropy field as 
the tunneling spectroscopy is obtained in the long-time 
limit. 

The tip and the substrate are assumed to be in local 
equilibrium which allows us to employ the fluctuation- 
dissipation theorem. For a general lesser/greater GF 
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FIG. 2: (Color online) (a) Typical dependence of the anisotropy 
energy Ft / s as a function of the high-energy cutoff D c . (b) Tran¬ 
sition energies £1 = £ 3/2 - £ 1/2 and £2 = £ 5/2 - £ 3/2 as a function 
of F for a S = 5/2 local moment with uniaxial anisotropy D = 0.7 
meV and £ = 0. Inset: The ratio between the two transition en¬ 
ergies. 


j7f <,/> / we can then make use of the relation J{ < ^ > (co) = 
( ±i)f(±a>)[-2lmjl r (a ;)], where F! is the corresponding 
retarded propagator. Although the effect of F c is negli¬ 
gible in the finite voltage bias regime, we write the full 
expression here for completeness, and also since Ft / s are 
easily obtained from these expressions in the limit V —> 0. 
We have 

r c = 2 ^-e- ileVt £ f [f(a>)lm[F; r (a>)F’ k (a> - eV)\ 

- (/(<u) - f(a) - eV))F; r (a>)lmF’ k (a> - eV)\dco. (9) 


Since the effect of F in the finite voltage bias regimes is 
discarded, we write D according to D = Do + Dy, where 
the subscript 0 (V) refers to zero (finite) voltage bias. The 
contribution from the tunneling current generates 

2T 2 r 

D c o=-^-lm'Y_ l J f(co)G’ p (co)G' k (co)da>-Fo, (10a) 

P k 

flf^)hn[G r p (co)G r k (co-eV)\ 

pk ' 

- (/(<£>) - f(a> - eV)^G p (oj)lmG k (oj - eVudco, 

(10b) 


where To is the zero voltage bias result. Here, we no¬ 
tice a few interesting issues with the current-induced 
anisotropy field D. In equilibrium, V = 0, the field Do 
is purely real which implies that it merely provides an 
overall energy shift of the spin spectrum. Moreover, in 
the nonsuperconducting case, the expression for Do re¬ 
produces fhe electronically mediated exchange interac¬ 
tions between spin moments at different points in space 
and/or time; e.g., see Refs. 25-29. Here, since we ad¬ 
dress a system with a single spin, this interaction can be 
understood as a self-interaction. 

We proceed with the calculations of the fields D and F 
by introducing the retarded GF for the surface electrons 
in the substrate which, in Nambu space, can be written 
as 


G!»=^ Gk 


p+r Q+r|(w) 
k u k 


CO + z6 + 


A 


(aj + id) 2 -E?\ A * a> + id-e k J' 


( 11 ) 


and analogously for the electrons in the tip. Sum¬ 
ming over the momentum, assuming energy indepen¬ 
dent density of electron states n sub , we then have 


G su b( a; ) “ W sub 


CO 


KubW - n sub 


-^w 2 -|A sub | 2 

Asub 
— [A sub | 2 


0(M |A sub |)A sub (<x>), 


Asub(w) 


(12a) 

(12b) 


and analogously for the tip GF. Here, we have intro¬ 
duced the notation A su b(w) = ln|l - 2o) sub /(w sub + C>c)l + 
msgn(a>)/2 where a> su b = y]co 2 - |A sub | 2 and D c is a high- 
energy cutoff, or half the width of the metallic band. 

In equilibrium, the Cooper pair correlations then give 
rise to the anisotropy from, e.g., the tip. 

In 

(13) 


1 - 2 - 


^sub 


Wsub + Dr 


dco. 


Fs=n 2 sub T: 


If 


/(w) 


|A sub | 2 sgn(a;) 
|A sub | 2 


CO- 


,2 


\ I though the integrand in this expression is logarithmi¬ 
cally small, it is non-negligible due to the integration 
over occupied states, which essentially spans the energy 
range [—D c ,0]. We can, moreover, see that this energy is 
finite, although decreasing with the high-energy cutoff, 
which is illustrated in Fig. 2(a). Hence, for a given band¬ 
width D c , the anisotropy constant oc (T s « sub |A sub |) 2 is es¬ 
sentially determined by the coupling strength T s , which 
varies exponentially with the distance between the local 
spin moment and the substrate, and analogously for the 
local moment and the tip. 

Here, we make contact with the experiments per¬ 
formed in Ref. 16, which involves a S = 5/2 local mo¬ 
ment, where the STM tip to molecular sample distance 
was varied in the tunnel junction. As an effect, the tran¬ 
sition energies drifted towards higher energies with de¬ 
creasing distance resulting from an increased uniaxial 
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anisotropy field acting on the local spin moment. By 
solving the effective spin Hamiltonian given by Eq. (7) 
as a function of for, e.g., a S = 5/2 local moment, we 
obtain the spectrum E n /2 = (n/2) 2 (D + !F), which gives 
the transition energies £i = £ 3/2 -£ 1/2 = 2 (D + 2 F) and 
£2 = E 5/2 “£ 3/2 = 4(D + 2F), as shown in Fig. 2(b), as 
well as a constant ratio £3 / £2 = 1/2; see inset of Fig. 
2(b). From this, we conclude that although the proxim¬ 
ity mechanism of the superconducting tip and substrate 
may not be the sole cause of an increased anisotropy with 
decreasing the tip-sample distance, the Cooper pair cor¬ 
relations within, e.g., the superconducting tip create a 
finite source for the anisotropy. As the effect from the 
Cooper pair correlations is essentially constant, which is 
crucial in the context, the increased anisotropy field with 
decreased tip-sample distance is provided by the expo¬ 
nential distance dependence in the coupling parameter 

bl¬ 
under nonequilibrium conditions, the modification of 
the spin excitation spectrum due to the tunneling current 


is not well understood at the moment since the tempo¬ 
rally fluctuating field induced by the Josephson current 
vanishes on average, but also since the field induced by 
the Josephson current is expected to decay with increas¬ 
ing voltage bias. Moreover, the field D c y only acts on 
S 2 and, hence, merely provides a rigid shift of the spec¬ 
trum. We, therefore, proceed with our discussion about 
the tunneling current and conductance. 


IV. THE TUNNELING CURRENT 
The tunneling current I(t) is derived from the relation 
I(t) = - edt^(dp a c pa ) = i2eIm^<Cp (T [cp (T ,'Wr]> / (14) 

p a pa 

with the electronic charge e and h = 1 , which leads to 




pk p'k' 

00 ' a"a’" 


r' (t)/ {t')T a " a "’ (t') c k'cr"' ^ 

(15) 

Mty P 'o"(ny eV{t - n }dt', 

(16) 


where T CTcr /(f) = Tq5 00 > + Ti a aa > ■ S„(f), whereas eV repre¬ 
sents the applied bias voltage by letting £ p(J —» e' po + eV/2 
and £^0 c' ki _ - eV/2. The two terms inside the expecta¬ 
tion value account for two different tunneling processes, 
of which the first gives the Josephson tunneling con¬ 
tribution and will be omitted throughout the rest of the 
paper since we are interested in the regime far from equi¬ 
librium where the Josephson correlations are negligible. 
The second term describes single-electron tunneling and 
may, analogously to previous studies on NM leads, be 
divided into three different parts proportional to Iq oc T 2 , 
h oc TqTi, and I 2 oc T 2 . 30-34 

Ip is the direct tip to substrate tunneling current. 


I 0 (t,V) =2eT 2 ReJ^£ [G^(f 

e mt- t ')d t ' / ( 17 ) 


pka 

■G> a (f',f)G< ff (f,f') 


moment, which can be seen from 


h(t,V) =2eT 0 T 1 Re^ f {a ao , ■ S(f) + <r 8B , ■ S(t'))e kv ^ 

, sj —OO 

pk 

00' 

x f G< a (t',t)G> ff ,(t,t') - G> a (t',f)G< a ,(t,f')]rff'. 

(18) 


This contribution is discarded, however, since it vanishes 
in the absence of spin-polarized currents. 30,31 . l 2 , finally, 
contains the spin-spin correlation function of the local 
magnetic moment coupled to the tunneling current. In 
general. 


hit 


,V)=2 e T 2 Re£ f 

1 *J — o 

pk 


G< a (t',f)G> a ,(f,f')v > (£l') 


and contains electron Green functions of the kind 
G p,a( f V) = ! '<Cpo( 0 c pff(O}'Ht ip / where *Htip determines the 
electron environment in the tip. Contributions from Iq 
are kept as a background current in our calculated re¬ 
sults. I\, on the other hand, does couple to the local spin 


‘ G P,a 


(t',t)G< a ,(t,t')x < (£t') 




(19) 


where *>(££) = (S(f)S(f')) and *<(££) = (S(t')S(t)). 
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The electron-spin projection onto the spin-correlation 
functions of the local magnetic moment amounts to 

g-qq' -<S(f)S(f ))-CTq'o = F {l-XaB +XgB + Xgg ) g ^ ^ ^ 

a/3 

( 20 ) 

where 

4>“ + ' + “ =<«|S z '-' + liS)<(i|S 2 ' + '-|n)P(£„)[l-P(E^)], (21) 

and the labels a, /S refer to the eigensystem {E a ,\a)} 
for the spin Hamiltonian in Eq. ( 6 ). In the follow¬ 
ing, we employ the decoupling ((dgd p )(t)(dj i d a )( f')) = 

,t), where g<(t,t') = i(dl(t')d a (t)) and 

g>(t,n = (- i)(d a (t)dUn > are the lesser/greater Green 
functions for the spin excitations. In the atomic limit, 

we have g^ > {t, t') = ( ±i)f(±E a )e~ tEa ^ t ~ t '\ where f(a>) is 
the Fermi function. 

As we consider stationary conditions in the setup, we 
Fourier transform our theory to the energy domain, us¬ 
ing, e.g., g(a>) = f git-t'y^-^dt'. 

A. Normal metal to superconductor junction 

For a normal-metal (NM) to superconductor (SC) junc¬ 
tion, the conditions differ for tip and substrate elec¬ 
trons/quasiparticles and, consequently, so does their re¬ 
spective Green functions. At the outset, considering the 
tip to be of a NM, the electronic structure is provided by 
Eq. (2a). The lesser and greater Green functions of the 
tip are then simply 

Gp'/'M =(±i)f(±co)6(£p - tv), ( 22 a) 

where f(co) is the Fermi-Dirac distribution function. 
Within the SC substrate, the quasiparticle structure given 
by Eq. (3) leads to lesser/greater Green functions, 

G^ > (o;) =(±z) ^\u k \ 2 f(±a>)5(E k + <n) + \v k \ 2 f (+a>)5(E k - «)}, 

(23a) 

where »k = Vl/ 2 (1 + e k/£k) and v k = ^1/2(1 - e k /E k ) 
are the coherence factors, and the quasiparticle energy 

E k = ^4 + IAsJ 2 . 

Following the approach in, e.g.. Ref. 35, we let 

L p -» fde ft ti p and L k -» dE k n sub E k / ^ - |A k | 2 in 

Eqs. (17) and (19), where fttip and n sub are the energy- 
independent tip and substrate electron/quasiparticle 
density coefficients, respectively, to obtain 

f(£ - eV) - /(£ + eV) 

I 0 (V) =4neT 2 0 n v n k dEE— - -- (24) 

d Ak V£ 2 -|Akl 2 


for direct tunneling and 

p- 2xle+X~J+X + a p 

I 2 (V)=4neT 2 1 n tip n sub £ j dEE ? ? 

Tp dA k y£ 2 -|A k | 2 

x {/(£- E a + E p - eV )-/(£-£„+ E p + eV) 

+ [/(£ + E a - E p - eV )-/(£ + £„- E p + eV) 

+ /(£ - E a +E a + eV) - /(£ -E a +E p - cV)]/(e)} 

(25) 

for spin-exchange tunneling. The remaining energy in¬ 
tegrals of the two expressions are solved numerically. 

Equation (24) provides a qualitative background 
shape for the current, as a function of bias voltage, that 
modifies under interaction with the local spin moment 
when (25) is added. Considering positive bias voltages, 
/(£ + eV) ~ 0 at low temperatures, while /(£ - eV) sud¬ 
denly jumps from 0 to 1 when £ and eV start to match. 
The tunneling current is, consequently, close to 0 while 
0 < eV < z\ k ; it quickly rises to a value determined by the 
quotient factor once eV ~ A k and transitions into a linear 
increase for eV > A k . In the conductance dI/dV(eV) spec¬ 
tra, the behavior is reflected in a sharp peak structure at 
eV = Ak, preceded by dl/dV ~ 0 A/V, and followed by a 
constant value. 

When the bias voltage is positive, the explicit expres¬ 
sion for h(V) suggests that a contributing tunneling 
channel opens up once eV > Ak + Ep — E a for transitions 
allowed by x z ap , X~ p , X„ p , i-e., that conserves angular 
momentum. In the dl/dV spectrum, this means that the 
initial peak at current onset is accompanied by smaller 
peaks for higher voltages corresponding to the spin ex¬ 
citation energies of the possible transitions. Apart from 
losing energy to the local spin by exciting it from the 
ground state, a tunneling electron may also gain energy 
from deexcitation of a thermally populated higher state. 
Such occurrences cause dl/dV peaks at lower voltages 
than the main peak. For negative voltages, the dl/dV 
spectrum is a mirror image with respect to eV = 0. 


B. Superconductor to superconductor junction 

Changing the STM tip from a NM to a SC, that is, 
using the quasiparticle structure given by Eq. (2b), we 
replace the tip GF by the ones given in Eq. (23) (replacing 
k p). Hence, using Eqs. (17) and (19), the direct and 
spin-exchange currents now become 
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(26) 


for direct tunneling, and 
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(27) 


for spin-exchange tunneling, where 0(x) is the Heavi¬ 
side step function, A p = A^ = A, and the superconducting 
phase difference <p = 0. The sign-alternating coefficients 
v+ = (-l)(« 2+J1 + 2 )/ 2 and v - = (—l)«(n+l )/2 c ]-i an g e s jg n for 
every other term starting with plus and minus, respec¬ 
tively. 

Despite the apparent added complexity, (27) behaves 
in much the same way as (25) with some qualitative dif¬ 
ferences. The onset of current by the applied bias volt¬ 
age no longer happens when eV ~ A, but instead occurs 
when eV ~ 2A since the step functions include an addi¬ 
tional pair potential. A, to the lower integration limit. 



E-Ef 


FIG. 3: (Color online) The bottom panel illustrates that the 
tunneling current onset happens at a bias voltage of V = A/e 
between a NM (orange) and a SC (turquoise). The DOS peak at 
onset carries over to the conduction spectrum. The top panel 
illustrates the same thing for a junction of SC:s where a bias 
voltage of twice the gap, A, is needed for conduction. 


The additional fractions and step functions in the math¬ 
ematical expressions for the currents also cause much 
sharper peaks in the dl/dV spectra in comparison with 
the NM tip setup. 


V. CONDUCTION SPECTRA AND ANALYSIS 


In contrast to a STM setup with normal-metal leads, 
the use of superconductors brings two main characteris¬ 
tic differences to the dl/dV spectra that we have touched 
upon. First, tunneling electron-induced spin excitations 
that are energetically within the superconducting gap 
of the system never occur until the bias voltage has 
passed the gap. Energy exchange between tunneling 
electrons and the local spin moment is, in other words, 
shifted to |V| = (A + Ep - E a )/e, for a NM-to-SC junction 
and | V| = (2A + Ep — E a )/e for a SC-to-SC junction, rather 
than | V| = (Ep - E a )/e, for a NM to NM junction. Though 
mathematically intrinsic, the physics picture to bear in 
mind is that tunneling from a SC at low temperatures 
only happens once enough energy is available to break 
up one Cooper pair. When tunneling to a SC, a single 
electron cannot occupy a SC low-lying state but must 
find a quasiparticle state higher in energy. The minimum 
energy cost for either event is A eV when one lead is a SC 
and 2A eV when two SC leads are used. Second, while 
inelastic scattering signatures in a NM-to-NM junction, 
of leads with flat density of states (DOS), appear as steps 
of increased conduction at the onset energies in the dl/dV 
spectra, *' 18,36,37 sCs produce peak structures followed by 
the usual stepped increase. These peak structures are left 
in the dl /dV curve as a trace by the underlying SC DOS, 
which exhibit pronounced coherence peaks at the end 
of the gap on both the occupied and unoccupied side. 38 
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FIG. 4: (Color online) (a) Schematic picture of the possible spin transitions for S = 1 with 0 or finite transverse anisotropy E. 
(b) Calculated conductance spectra of a SC-SC junction for local spin moment S = 1 at T = 1.2 K with parameters D = 1 meV, 
Af ip/sub = 0-5 meV, 7\ = 0.3 Tg, and varying E. (c) Same as in (b) for a NM-SC junction at T = 0.5K. 


Just as the bias voltage provides enough energy for an 
additional conduction channel to open, either the occu¬ 
pied or unoccupied states are inevitably at peak density. 
The conduction is momentarily high and falls off once 
the bias has passed the peak. See Fig. 3 for an illustrative 
description of these tunneling properties. 


A. Spin 1 magnetic molecule 

Three quantum states \m z = —1,0,1) exist for a local 
magnetic moment of S = 1 and a finite axial anisotropy 
D generates two energy eigenvalues to the spin Hamil¬ 
tonian 'He, if the transverse field E and the external mag¬ 
netic field B are absent. A positive anisotropy D > 0 
meV will cause the eigenstates \m z = ±1) to lie E± = D 
meV above the state \m z = 0), thus favoring a low-spin 
ground state. A negative anisotropy D < 0 meV will, 
instead, favor the high-spin state since E± < Eg. 

The left panel in Fig. 4(a) schematically shows the 
possible spin transitions in the case D > 0. Hence, feed¬ 
ing the energy corresponding to |D[, in addition to the 
energy needed to overcome the superconducting gap(s), 
into the system allows the local spin moment to undergo 
transitions between its ground and excited states, which 
is clearly illustrated in the bottom traces of Figs. 4(b) and 
4(c). Here, additional conductance channels emerge at 
eV = ±2| A| + D in the SC-SC junction and eV = ±| A| + D in 
the NM-SC junction, respectively, caused by the inelastic 
scattering. In this case, we have used the values A = 0.5 
meV and D = 1 meV. Qualitatively, the conductance spec¬ 
trum in Fig. 4(b) agrees well with the experimentally 
obtained conductance in Ref. 13. 

The dl/dV curves of the SC-SC and the NM-SC junc¬ 
tions in Figs. 4(b) and 4(c) are in stark contrast to each 
other in terms of peak width. The SC-SC peaks are very 
sharp even though the calculations were done at a tem¬ 
perature of 1.2 K as opposed to 0.5 K for the NM-SC 
case. The difference is to be expected to some extent 
since two DOS coherence peaks match up at the onset of 
any new conduction channel to give a very conductive 


SC-SC junction for a narrow voltage span. In contrast, 
the conduction for a NM-SC junction, where one lead 
has a flat DOS, differs less at onset voltage in compari¬ 
son to higher voltages. While this reasoning will explain 
a noticeable difference, the huge discrepancy found in 
our calculations indicates a failure of theory to handle 
peak widths in the SC-SC situation. 

For E + 0, the eigensystem of the local spin is mod¬ 
ified to E ± = D ± E, |E ±1 > = [\m z = -1) ± \m z = 1)]/ V2, 
which breaks the degeneracy and separates |E + i) from 
|£_i) by 2E in energy. The spin changing transitions, 
e.g., (£o|S + |£_) and (£o|S_|£ + ), therefore occur at dif¬ 
ferent energies, as illustrated in the right panel of Fig. 
4(a), and we expect conductance signatures at the volt¬ 
age biases \eV\ = A su t, + D ± £ for the NM-SC setup and 
at |eV| = A t ip + A su b + D±E in the SC-SC case, which is 
readily seen in Figs. 4(b) and 4(c), respectively. The 
dashed and dotted lines trace the actual progression of 
eigenvalue differences with respect to increasing val¬ 
ues of the transverse anisotropy. In addition, because 
the Fock states \m z = +1) are coupled, the tunneling 
current also facilitates spin-preserving transitions be¬ 
tween the states |£ + i) and |E_i). Inelastic signatures 
between these higher-energy spin states are expected 
to appear on both sides of the main coherence peaks. 
At \eV\ = Atip + A su b + £, the dotted line leaning towards 
the right in Fig. 4(b) traces peaks from excitations 
|£_) —> |E+). The dotted line leaning towards the left 
traces the barely visible in-gap peaks, indicated by an 
arrow in the middle curve, at \eV\ = A t ; p + A su b - £ from 
deexcitations |E+) —» |£_) that assist electrons in tunnel¬ 
ing. The higher-energy states of the local spin reveal 
themselves in this manner since they are thermally pop¬ 
ulated enough at kgT oc O.lmeV to support transitions. 
The NM-SC dl/dV curves of Fig. 4(c) are calculated 
at a lower temperature that populates the higher-spin 
states less, which in turn prevents a clear signature from 
|E + i) fc; |£_i) transitions. 

The apparent difference in amplitude between the 
transitions |£o) —> |E±i) and |E+i) —» |E+i), which is leg¬ 
ible from Fig. 4(b), can be understood in terms of the 
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population factors P a p. For D > 0 and small E ± 0, the 
populations P± of the states |£+i) are both close to 0, 
such that, e.g., (1 — P + )P- becomes small. The popula¬ 
tion Pq for the state |£q) is, on the other hand, close to 1 

1.5 

d 

> 1 
±5 

(a) 

A 


has-o-'e 


A 

which leads to relatively large products (1 - P±)Po- Note 
also that as P_ gets larger for greater values of £, while 

P + gets smaller, the (1 -P+)P- peak gets bigger. At the 

0.5 



fW =0.065 






=0.026 




gets slightly bigger, even though the low initial value of 
P + prevents any considerable changes. 


B. Spin 5/2 magnetic molecule 

Next, we turn our attention to the spin S = 5/2 
system in order to connect to recent experimental 
observations. 1 * 1 For £ = 0, the eigensystem consists of 
the doubly degenerate states \m z = ±m/2), m = 1,3,5, at 
energies E ±m / 2 = Dm 2 /4, and, with a positive (negative) 
uniaxial anisotropy, D > 0 (D < 0), the system acquires a 
minimal (maximal) spin state | ± 1/2) (| ± 5/2)). 

In Fig. 5(a), we plot the calculated SC-SC junction con¬ 
ductance for varied populations of the states \m z = ±3/2) 
in the absence of transverse anisotropy, £ = 0. We infer 
that our model calculations reproduce the experimental 
observations with excellent agreement. Here, we as¬ 
sume that the pairing potentials of the tip and substrate 
are equal, A t i p / su b = A ~ 1.35 meV, neglect possible super¬ 
conducting phase differences, and use a positive uniax¬ 
ial anisotropy D = 0.7 meV. Analogously to the previous 
case, the conductances display strong coherence peaks 
at eV = ±2A, which are perfectly replicated at the volt¬ 
age biases \eV\ = 2A + 2D for the inelastic spin transition 
I m z = ±1/2) -»| m z = ±3/2). 

We, furthermore, notice the conductance peak emerg¬ 
ing at voltage biases \eV\ = 2A + 4D for an increased pop¬ 
ulation of the first excited states \m z = ±3/2). The con¬ 
ductance peak is a signature of the inelastic transition 
| m z = ±3/2) —» | m z = ±5/2) and its characteristics can be 
quantified by using the expressions in Eqs. (20) and 
(21). As the matrix elements for raising and lowering be¬ 
tween the states | m z = ±3/2) and \m z = ±5/2) are always 
finite in the present setup, the emergence of the conduc¬ 
tance peak strongly depends on the population of these 
states. When the ground state is heavily populated, both 
| m z = ±3/2) and \m z = ±5/2) are largely unpopulated and 

the factors P. 3.5 are vanishing. This scenario remains 

±2 = 2 

valid for small charge currents through the system, as 
well. For increasing charge currents, however, popu¬ 
lation density is expected to accumulate in the states 
\m z = ±3/2) as they are excited with a faster rate than their 
corresponding decoherence times. Accordingly, upon 
populating these states, the factors P + 3 ± 5 become finite 
which leads to the transitions \m z = ±3/2) —» \m z = ±5/2) 
contributing additional channels for conduction. In this 
fashion, we reproduce the effect of pumping which is 
obtained experimentally by decreasing the distance be¬ 
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FIG. 5: (Color online) (a)-(c) Calculated SC-SC conductances 
for a spin S = 5/2 system under (a) varying population of the 
states | m z = ±3/2) for D = 0.7 meV, E = 0, (b) varying uniaxial 
anisotropy D, for E = 0, and (c) varying transverse anisotropy 
E, where D = 0.7 meV and P ± 3/2 = 0.16. Other parameters are 
Atip/sub = 1-35 meV, T = 1.2 K , 16 and Tj = 0.3T 0 . (d)-(f) Cor¬ 
responding conductance spectra for a NM-SC junction for pa¬ 
rameter values T = 1 K and T\ = 0.1 To, while other parameters 
are as in (a)-(c). 


tween the scanning tip and the sample. 

Figure 5(d) illustrates the corresponding conductance 
spectra for a spin-5/2 magnetic molecule trapped within 
the gap of a NM-SC junction. Once again, no qualita¬ 
tive differences are obvious from the SC-SC case, ex¬ 
cept for wider peaks and earlier onset, at bias volt¬ 
ages \eV\ = A for the main conductance peak and at 
\eV\ = A ± 2D for the \m z = ±1/2) —> \m z = ±3/2) transition. 
With higher population numbers of the state \m z = ±3/2), 
which are motivated if the local spin mainly dispenses 
excitation energy and angular momentum to the SC 
substrate through the relatively slow spin-phonon cou¬ 
pling to allow for pumping, signatures from the inelastic 
| m z = ±3/2) —»| m z = ±5/2) transition are revealed. 39 

The plots in Fig. 5(b) and 5(e) show the evolution of 
the IETS spectra as a function of the effective uniaxial 
anisotropy D, which may be thought of as the sum of 
the intrinsic molecular anisotropy and the Cooper-pair- 
induced anisotropy; cf. the model in Eq. (7). The shift 
to higher energies of the inelastic peaks is expected from 
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the previous discussion in Sec. III. In the experiment, 
the STM tip was brought closer to the sample which is 
expected to generate an exponential growth of the tun¬ 
neling current since the coupling between tip and sam¬ 
ple varies exponentially with distance. Here, as we do 
not attempt to model the whole experimental setup, but 
rather investigate the effects of changes in the anisotropy, 
we have not included this exponential variation of the 
anisotropy. In the case of a NM-SC junction, we may 
think of the tip as superconducting while the substrate 
is normal metallic for a feasible physical setup where an 
increase of the effective D follows when the STM tip is 
brought closer to the sample. 

For a finite transverse anisotropy, E 0, a peak can be 
seen to rise along the dash-dotted line in the SC-SC panel 
of Fig. 5(c) as the value of E gets bigger. To explain the 
appearance of this peak, we look at how the spin states 
modify simultaneously to form linear combinations of 

the kind \E ±m ) = X!,„=i, 3,5 / 2 l m - = ±n / 2 )- The six spin 

states are still doubly degenerate on three energy lev¬ 
els, but there is now a finite probability that a transition 
from the lowest state, e.g., |Ei) weighted on \m~ = 1/2), 
to the highest, e.g., IE 3 ) weighted on \m z = 5/2), occurs 
despite seemingly violating conservation of angular mo¬ 
mentum. Consequently, increased values of E distribute 
density among the Fock states to allow for transitions 
with A m z = +1 between any of the available states. A 
schematic picture of the added transition possibilities 
for nonzero transverse anisotropy £ is given in Fig. 7. 
With different values of E, the spin state energy levels 
also shift relative to each other, which is reflected in the 
peak positions of Fig. 5(b). For example, at just over 
E ~ 0.2 meV, e.g., the |£i) —» IE 2 ) and IE 2 ) —» IE 3 ) transi¬ 
tions clearly cross in energy. 

The characteristics of the SC-SC conductance spectra 
translate, once again, to the NM-SC case for finite values 
of £ since both systems share the local spin structure; 
see Fig. 5(f). Spectral details of the internal workings 
are, however, easily lost in the thermal broadening of 
the transition peaks. 

C. Influence of the temperature 

We notice in Eq. (27) that there is an increased degree 
of detail in the tunneling current and, hence, the conduc¬ 
tance spectra, at elevated temperatures, provided that 
we remain below the critical temperature. In Figs. 6 (c) 
and 6 (d), these details are reflected in the conductance 
spectra as additional peaks that rise with increased pop¬ 
ulation of states higher in energy. For the in-gap peaks, 
there are two mechanisms responsible, which are of sim¬ 
ilar origin. 

To begin with, consider the S = 5/2 system, whose 
spectrum for finite D > 0 is provided schematically in 
Fig. 6 (a). For low temperatures and small currents, the 
excited states are expected to be more or less unoccupied, 
while the ground state, covered by a finite portion of the 



FIG. 6: (Color online) (a) Schematic spectrum of a spin S = 5/2 
SC-SC system in the atomic limit subject to different condi¬ 
tions. The blue, red, and green traces in the right panel rep¬ 
resent the Fermi function at temperatures T = 1.2/4/6 K. (b) 
Schematic DOS illustration depicting temperature effects for a 
S = 1 system, (c), (d) Conductance spectra for the SC-SC setup 
at temperatures corresponding to those given in (a) and (b) for 
a S = 1 and S = 5/2 system, respectively. Other parameters are 
as in Fig. 5. 


Fermi function (blue), is occupied. Under these con¬ 
ditions, we retain the previous conductance spectrum, 
which is reproduced in the bottom trace of Fig. 6 (d) for 
convenience. 

At elevated temperatures, the thermal excitation en¬ 
ergy is sufficient for the excited states to be partially 
occupied; see Fig. 6 (a), red and green Fermi functions. 
The local spin can then undergo transitions not only from 
lower to higher excitations, but also from higher to lower. 
In the former case, the spin moment has to absorb energy 
from the tunneling current, hence, the voltage bias has 
to be sufficiently large to assist such a transition, e.g., 
\eV\ < 2|A| +£ 3 / 2 -£ 1 / 2 , where £ 3 / 2 -£ 1/2 >0. In the latter 
case, however, the spin moment is already thermally ex¬ 
cited and can undergo deexcitation processes at energies, 
e.g., \eV\ < 2|A| +E 1/2 - £ 3 / 2 , where £ 1/2 - E 3/2 < 0. Then, 
the spin moment emits the energy quanta E 3/2 - E 1/2 into 
the tunneling current, a process that opens a new chan¬ 
nel for conduction which is expected to be seen within 
the gap of the conductance spectrum. This is indeed 
the case, evident in the conductance traces calculated for 
T = 4 and T = 6 K in Fig. 6 (d). The peak near equi¬ 
librium corresponds to the transition | ±3/2)(±5/2| since 
4D = 2.8 meV — 2|A| = 2.7 meV. Given that the thermal ex¬ 
citation energy kgT is greater than the difference between 
the two energies, ~ 0.3 and 0.5 for T = 4 and 6 K, respec¬ 
tively, there is room for these deexcitation processes near 
equilibrium. The second in-gap feature corresponds to 
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FIG. 7: (Color online) (a) Spectrum of a spin S = 5/2 SC-SC system in the atomic limit subject to different conditions, parametrized 
by the anisotropies D, E and external magnetic field B = Bz. (b) Corresponding conductance spectra for a NM-SC setup. Unspecified 
parameters are as in Fig. 5. 


the emission resonance for the transition | ± l/2)(±3/2|, 
emerging equidistantly from the superconducting coher¬ 
ence peak as its corresponding absorption resonance. 

While spin deexcitation assisted tunneling certainly 
accounts for some of the in-gap features, we continue our 
discussion by looking at the spectra for S=1 given in Fig. 
6(c). Only one excitation peak exists alongside its deex¬ 
citation signature close to the main coherence peak. Yet, 
there is still an emerging three-peak structure forming 
with higher temperatures in the center of the supercon¬ 
ducting gap. A second mechanism is clearly at play and 
it can be explained by looking at the schematically drawn 
DOS in Fig. 6(b). For low temperatures, the Fermi func¬ 
tion (blue) occupies the subgap states only, but for higher 
temperatures, the Fermi function (green) stretches all the 
way to the overgap states that become slightly filled. The 
conductance is consequently nonzero for low-bias volt¬ 
ages through direct spin-preserving tunneling. Once 
the voltage bias passes the E+i - Eo difference, spin-flip 
tunneling, which excites the local spin, may also occur, 
adding a peak at eV = D meV. In-gap resonances were 
recorded for a S = 1 Mn-phthalocyanine at T = 4.5 K, 13 
and we believe that the in-gap resonances considered 
here can, at least partly, explain these observations. 


D. Spin 5/2 magnetic molecule under external magnetic 
field 

In order to explore additional aspects of the conduc¬ 
tion spectra for the S = 5/2 system, an external magnetic 
field is introduced to break up the twofold degeneracies 


that the anisotropy fields D and E are unable to. Figure 
7(a) pictures the expected behavior for three different 
magnetic field intensities in the z- direction. For E = 0, 
two smaller peaks emerge around the main coherence 
signature, at eV = 2A, equidistant on both sides as a re¬ 
sult of inelastic emission and absorption between the 
Zeeman split ground states, \m- = -1/2) ±5 \m z = 1/2). 
Forking off as the magnetic field increases at about 
V = 4 mV are two peaks that signal transitions between 
<: m z = 3/2|S+|m z = 1/2) and (m z = -3/2|S_|m z = -1/2). 
These transitions differ in energy because the pair of 
ground states are Zeeman split by a different amount 
than the first excitation pair of states. Note that the tran¬ 
sitions occur between uncoupled basis states since £ = 0, 
which limits the number of possible excitation paths to 5, 
as schematically illustrated in Fig. 7. Indications of tran¬ 
sitions between the first and second pair of excitation 
states are absent unless the effects of pumping are repli¬ 
cated as done in the previous example with no external 
magnetic fields. 

For £ > 0, the basis states once more couple to form 
eigenstates to the spin Flamiltonian. This is reflected 
in the dl/dV plots of Fig. 7(a) for increased magnetic 
fields in the z-direction as a branch off of the transition 
signature, at V ~ 4 mV, into four peaks rather than the 
previous two for B z = 0. Any transition |E ± i) —» jE+ 2 ) is 
hence sufficiently probable to yield a visible peak in the 
conductance spectra. For E ~ 0.2 meV, we even begin to 
see four distinct peaks split off, in step with the magnetic 
field, that originate from excitations between the ground 
and the second excited states at V ~ 7 mV. In theory, 
transitions are now allowed between all spin states of 
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the magnetic molecule, once D > 0, E > 0, and B z > 0, 
even though thermal populations for all but the two 
lowest-energy states are so small that excitations from 
higher states are rare occurrences; see upper left corner 
of Fig. 7 for a diagram of the possible excitations. 


In Fig. 7(b) we look at the system under equal cir¬ 
cumstances regarding the external magnetic field for 
the NM-SC setup. Unfortunately, reasonable magnetic 
fields separate the peaks from different spin transitions 
less than the thermal width which somewhat obscures 
details. Features of the underlying peak structure can 
still be made out as additional humps form with stronger 
magnetic fields, but aside from resolving in energy, it is 
possible to draw conclusions based on amplitude. For 
E = OmeV, a single peak appears to form at just over 
V ~ 2 mV with an amplitude that is strongly depen¬ 
dent on the magnetic field. What appears to be one 
signature is really two peaks that separate for stronger 
fields. The peak moving towards the right, originating 
from |nu = 1/2) —> |3/2), quickly dies off as \m z = 1/2) 
becomes less populated at the low temperature, while 
the peak moving towards the left, from the excitation 
| m z = -1/2) —> | - 3/2), gains amplitude as \m z = -1/2) 
becomes more populated. In this way, an external mag¬ 
netic field can assist to increase amplitude for some tran¬ 
sitions. 


A magnetic field in the x- and y- directions similarly 
splits up the degenerate energy levels of the local spin. 
Around the main coherence peak, signatures from both 
absorption and emission can be seen when the spin leaps 
in energy between the separated ground states. Starting 
at approximately V ~ 4mV, we see in Fig. 8 how the 
excitations |Ei(-i)) —> |£2(-2))/ which share energy, pro¬ 
duce a peak that divides into four when all transitions 
|£y_i) —> |£ 2 ,— 2 )/ at different energies, are allowed with 
the magnetic field, even though £ = OmeV. In compari¬ 
son with the B- 4 0 setup, these peaks break apart along a 
bent path rather than following a straight line. The most 
frequent transitions are also those with higher energy as 
opposed to those with lower energy. When the trans¬ 
verse anisotropy E is turned on, the conductance spectra 
look quite different as two peaks do not seem to sepa¬ 
rate, while the other two go off in opposing directions to 
effectively form a structure of three peaks. 


Providing an external magnetic field adds a compli¬ 
cation to the measurements, since the superconductivity 
in both the substrate and tip becomes quenched under 
too strong fields. This problem can, however, be over¬ 
come by changing to a tip/substrate material that is less 
sensitive to magnetic fields, e.g., NbTi, Nb 3 (Sn,Ge,Al), 
and MgE> 2 , 111 12 which are known to maintain their su¬ 
perconducting phase for fields as strong as 10-30 T. Our 
predictions made for fields up to a few T are therefore 
safely within the realm of feasibility. 
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FIG. 8: (Color online) dl/dV spectra of a S = 5/2 system under 
external magnetic fields in the x- direction with increasing field 
density. Only the transverse anisotropy E differs between the 
two diagrams. 


VI. SUMMARY AND CONCLUSIONS 

We argue that our simple model of a superconducting 
STM, holding a paramagnetic molecule within its gap, 
generates a differential conduction spectra that matches 
up very well to experimental data, taken of, e.g., Fe-OEP- 
C1 and Mn-phthalocyanine. The model notably captures 
peak signatures in the tunneling conductance from inter¬ 
actions with the local spin that reference to the sum of the 
tip and substrate pairing potentials rather than zero-bias 
voltage. We are also able to mimic the effects of electron 
pumping by introducing a uniform potential shift such 
that the excited spin states thermally populate to reveal 
peak imprints of transitions among them. The success 
up to these points leads us to infer that the key mech¬ 
anism behind the experimental conductance features is 
exchange interaction between tunneling electrons and 
the local spin moment. 

Our model does not include direct exchange between 
the local spin and the superconducting substrate, which 
will generate states within the superconducting gap, 
since we argue that the separating ligand cage weakens 
this interaction such that, e.g., Shiba states move close 
to the dominating coherence peak. Our model does, 
however, capture the effect of exchange between the lo¬ 
cal spin moment and the Cooper pair correlations that 
generates a finite contribution to the uniaxial anisotropy 
that acts on the local spin moment and increases with de¬ 
creasing distance between the superconducting tip and 
the sample. This effect, therefore, offers an explanation 
for the increased anisotropy observed in experiments. 

With the freedom to explore parameter space, we con¬ 
sider different magneto crystalline anisotropy values as 
well as the effects of an external magnetic field. For 
S =5/2, the axial anisotropy field directly determines 
level spacing between spin states, while the transverse 
anisotropy field, apart from slightly shifting the energy 
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levels, couples the spin basis states to allow for tran¬ 
sitions which are otherwise prohibited by conservation 
of angular momentum. An external magnetic field re¬ 
moves spin state degeneracies and provides a rich con¬ 
ductance spectrum. 

We have also considered temperature effects up to the 
critical temperature and shown that both direct thermal 
excitations of the local spin, as well as thermal popu¬ 
lation of states above the superconducting gap, cause 
in-gap peaks in the conductance. Deexcitation of ther¬ 
mally populated higher-spin states assist electron tun¬ 
neling at voltage biases lower than the superconducting 
gap, while thermally occupied states above the super¬ 
conducting gap give a nonzero conductance contribu¬ 
tion that peaks at zero voltage bias and when it matches 
up with the excitation energy of the local spin. These 
mechanisms may partially explain the observed in-gap 
resonance of Mn-phtalocyanine. 

An extended experimental study of the system could 


benefit from the use of an external magnetic field. The 
main argument for long spin excitation lifetimes is that 
deexcitations with an energy release between 0 < e < 2A 
fail to split up Cooper pairs and facilitate particle-hole 
creation. A magnetic field immediately produces a large 
peak that separates from the main coherence peak due 
to transitions between the no-longer degenerate ground 
states. The energy of this excitation varies with the 
strength of the magnetic field starting from OmeV and 
upwards. 
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